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1 

Descrintinn 

METHOD AND APPARATUS FOR SYNTHESIZING LEADS 
OF AN ELECTROCARDIOGRAM 

5 Techni^l Fi,!,) 

This invention relates generally to a sysrem for leads „ f 

an etecttocardrogram. and more specifically, ,o a method and an appLTfor 
synthestzutg leads based on developing a pauen-specific transform. 

10 BackrvoMnri »f a. rrrr-mim 

• • ^ '"^o^ogram (ECG) is an imponanr tool for monitoring 
hean aenvtrv and diagnosing heart condidons. Tne ECG is a recording ofZt 

^^T 1 ^'""- ^^^^cansearhehearrto™ 
m °^ bl °°< » * TO mroughont the bodv. This 
elecmcai aomty is spontaneously generated. As the cells within a. hean,'chan g e 

Sl'dT™ '° 1 P ° SW ™ POtemial (^^onX the muscles 

wtthtn the bean contract. Convetaely, „„„ ^ ^ , 

^^ T Tnrr^d• (repo ' aIization) • *• — - remra ^ 

state. The penodrc contracnon of the bean causes the pumping action. This 

rr:^ « *» — ™ £U 

a pauent s ECft vanous cardiac abnonnalides. such as ischemia, can be detected 
nlaced ™ * w elM °T i ° £ «" hean <="> be monitored by elecn-odes 

eLttT "Jf" ° f * "'P 01 ^ - "Polarize, the 

25 ST? ^ 0tt ^ ^ »' ^ ^ varies. Each conn4on o the 
j£ 0 ~ B b - ) a " TeSP °° ds t0 « «-*« depolarizaoon/repoiarizauon 

ootenrW m h"" Sta '" illrd """""^ A lead is the electrical 

between «, points on me body surface or (2) between one point and 

1Z h "T"* ^ 6 Sh0W 1 '^cardiogram 

d*playmg the twelve 1 ends. The twelve leads are referred to as L H. EL aVR. 

r s i T' 2 ' ' V4, V5 ' Md Vfr tet 561 ° f «* leads are known as 

^ - «* >* - shown in Figure 7. Tie fourth hntb etaode. P Ld on me 
fl* J» . a common ground for me enure system. Lead I is the potential 
duTerence between the left arm and the right arm. Lead H is the potential 
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difference between the left leg and the right arm. Lead IH is the potential 
difference between the left arm and the left leg. The other three limb leads - 
augmented voltage right arm (aVR), augmented voltage left arm (aVL), and 
augmented voltage left foot (aVF) - are the potential difference between their 
5 respective limb and the average potential of the other two limbs. Figure 7 
illustrates the measurement of aVR. The six limb leads are related 
mathematically. If any two of leads L II, or IH are given, then the other four limb 
leads can be calculated. 

-_ Leads V lf V 2 , V 3 , V 4 , V5, and Vg are known as chest leads. The 
10 electrodes for measuring these leads are positioned as shown in Figure 8. These 
leads measure the potential difference between the electrodes and a common 
reference known as the Wilson Central Terminal. The Wilson Central Terminal is 
formed by connecting the electrodes on the right arm, left arm, and left leg 
through a resistor network to a common point 
15 To aid analysis, the ECG complex, which is the portion of the ECG 

associated with the electrical activity, is divided into three components: P wave, 
QRS complex, the T wave. The P wave corresponds to the depolarization of the 
atria, the QRS complex corresponds to the depolarization of the ventricles, and 
the T wave represents the repolarization of the ventricles. The repolarization of 
20 the atria occurs simultaneously with the depolarization of the ventricles, but is 
typically masked by the strong electrical signal generated by the ventricle 
depolarization. 

In order to generate a twelve lead ECG, prior an techniques placed 
10 electrodes on the patient as shown in Figures 7 and 8. However, placing so 

25 many electrodes is very time consuming. Also, electrodes have a tendency to 
move on the body's surface or to fall off the surface. Such electrode movement 
typically accompanies patient movement The electrode movement results in 
inconsistent ECG signals, which make patient diagnosis difficult The use of ten 
electrodes and the ten wires is uncomfortable for a patient 

50 Consequently, only two electrodes are typically placed on a patient 

Although the use of two electrodes is sufficient to monitor arrhythmia, it is 
■ insufficient to diagnosis problems such as ischemia. When an episode of silent 
ischemia occurs in a patient then addition electrodes are attached so that the 
twelve lead ECG can be generated. However, because of the time needed to 

!5 attach the additional electrodes, the episode of ischemia may be over by the time 
the twelve lead ECG is generated. 
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A further problem in generating the twelve lead ECG is that the 
placement of the electrodes at the standard positions may be difficult, if not 
impossible, due to a patient injury or suture. Thus, it is impossible to generate the 
standard twelve lead ECG for these patients. 

Prior systems have attempted to synthesize the twelve lead ECG 
using fewer than the standard number of ten electrodes. However, all prior 
systems have produced unacceptable results. Some prior systems have attempted 
to synthesize the twelve lead ECG based on data gathered from orthogonal leads 
An orthogonal lead system measures the electrical activity of the heart in an XYZ 
coordinate system These prior orthogonal systems have produced unacceptable 
symheazed leads. Moreover, the use of orthogonal leads implies that electrodes 
are placed on a patienrs chest and back. Such a placement of leads is 
uncomfortable for the patient Some prior systems have attempted to provide a 
population-based synthesis. The population-based systems categorize patients 
based on several factors such as age, sex, body build, or electrical orientation of 
the heart The systems then synthesize certain leads based on data gathered from 
other leads and patient category. These systems, however, have produced 
unacceptable results. 

It is desirable to have a system that would produce an acceptable 
^thesis of the twelve lead ECG data gathered from less than ten electrodes. 
Furthermore, it is desirable to have a system that would produce an acceptable 
synthesis of a lead based on collecting of minimal amounts of actual lead data. It 
is desirable that such a system iiiaximizes patient comfort 
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25 Summary gfj^iffiaajfln 

It is an object of the present invention to provide a method and 
apparatus for synthesizing an ECG lead based on the actual data received for 
otner leads. 

It is another object of the present invention to provide a method and 
apparatus or generating a standard twelve lead ECG based on receiving data 
tram three leads. 

It is another object of the present invention to provide a method of 
segmenting the ECG complex to improve the methods of synthesizing lead data. 

These and other objects, which will become apparent as the 
invennon is more fully described below, are obtained by providing a method and 
apparatus for synthesizing lead data. In a preferred embodiment of the present 
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invention, the method for synthesizing data for a given lead comprises the steps of 
selecting a plurality of base leads, gathering a first set of ECG data from the 
patient for the base leads and for the given lead for an interval corresponding to at 
least one ECG complex, generating a transformation based on the first set of ECG 
5 data, gathering a second set of ECG data from the patient for the base leads for 
an interval corresponding to at least one ECG complex, and applying the 
transformation to the second set of ECG data to effect the synthesis of the data 
for the given lead. 

Brief Description of the Drawi^ 

10 Figure 1 is a block diagram showing Phase 1 and Phase 2 of a 

preferred embodiment of the present invention. 

Figure 2 is a block diagram showing the components of Phase 1 in a 
preferred embodiment. 

. Figure 3 is a block diagram showing the components of Phase 2 in a 
15 preferred embodiment 

Figure 4 shows the ECG complex divided into three segments. 
Figure 5 shows the ECG complex divided into eight segments. 
Figure 6 shows a sample twelve lead ECG. 

Figure 7 shows the placement of the electrodes used to generate the. 

20 six limb leads. 

Figure 8 shows the placement of the electrodes used to generate the 
six chest leads. 

Figure 9 is a flow diagram of calculate coefficients routine. 
Figures 10A through 101 comprise a flow diagram for the regression 

25 routine. 

Figure 11 is a flow diagram of the synthesized lead routine. 

Detailed Desq^nn Q f the Inven^ 

In a preferred embodiment of the present invention, a lead is 

30 synthesized based on a patient-specific transformation of data acquired on three 
base leads. The first phase of the present invention involves the generation of the 
transformation data. The transformation data is generated based on data from 
the three base leads and data from the lead to be synthesized. In essence, the first 
phase is a learning process in which a' system in accordance with the present 

J5 invention learns" the relationship between the base leads and the lead to be 
synthesized for a specific patient Once this relationship is "learned," then 
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transformation data for .he patiem is generarad. In me second phase, ,he present 
~ data from d,. three base leads from rhe patient. It the 

l,h hT 00 " 10 * *■ 10 8eneme S ~ ed *» f °' *• l=ad .He 
symheszed data a used to antuyze .he electrocardiographic activity. 

lea* In PrCf ! Ired embodiment - "* P"*«» "»™»on oses standard 

eads L n. and V, as d.e toee base leads. However, other standard leads anJ 
non-srandard leads wffl produce accept results, m general, .he greararl 
^dependence of the mree base leads, the better me^heTI £ 

10 I^ h W0Uld ^ ,0 ~= - dTe 

° nhOSOnal P "~ ~* 41 " placemen, o. 

In a preferred embodiment me presen. invention divides Ibe ECG 

Z el , Se r ntS - * ^ i,,Venti0, ■ 8 ™ ratts ^"nation dara for 
each segment When synthesizing the !ead dara, the present invention determines 
15 m whtch segment each dara poin, is in and applies the appropriate e^em 
transformation to the base leads. segment 

HB" 1 is » block diagram showing Phase 1 and Phase 2 of a 
preferred embodiment of me syndesis sysram. In Phase 1 (110) 2 JLZ 

20 hnear transformation. The inputs . Phase 1 are me mgitized'dau gen raL 
^.tite three b^e leads (Lead A, Lead B, and Lead C) and me dltTme 
bad » be synmea^d (Lead YL Phase 1 transfers me segment definitions £ 
rhe coefficemt , » Phase 2 (120). Phase 2 generates synthesized dam for Lead Y 
Tneoutp„ t ofPhase21sLeadY-. As Phase 2 receives data from me ba^Teal' 
25 it generates the dara for Lead Y'. •> me oase leads, 

n™*— f!T 2 15 a bIock dia 8 ranl 1» component of Phase 1 in a 

preferred embodunent He components illustrated are signal average (210) 

30 a^Ko^! , , T OTOTSe «™ - 

average ECG complex for each of tite input leads. In signal analysis, me use of an 
average stgnal rands to minimize the effect „, noise oaZ signTln a p^ 
embodmient, ran R-R cycles are inpu. for each ,ead. AltcZ^. TZt 
need on*- consider me ECG complex itself and no. me entire R-R JT 

* E^ZT " ** " ^ «~ and secern 

ECG complex component In an alramara embodiment oiher signal anaL 
rachmques may be used, such as, using me median point of each of me ranges 
as a subsntura for me average value. In me .amove pacemalter spike ^otnt 
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the system determines if there is a pacemaker spike present in the averaged 
signals. The techniques for determining whether a pacemaker spike is present are 
well-known, such as detecting a high-frequency waveform corresponding to a 
ventricular or atrial spike. The outputs of the remove pacemaker spike 
5 component are the averaged signals with any spike removed. 

In use segment ECG complex component (212), the system 
generates a segmentation of the ECG complex. In a preferred embodiment, the 
segment ECG complex component divides the ECG complex into three segments 
delimited by [1] the beginning of the ECG complex, [2] the beginning of the QRS 

10 complex, [3] the end of the QRS complex, and [4] the end of the ECG complex. 
Figure 3 illustrates the three segments of the ECG complex. In an alternate 
preferred embodiment, the segment ECG complex component further subdivides 
the segments of the ECG complex. As illustrated in Figure 5, the ECG complex is 
divided into eight segments. The eight segments are divided by [1] the start of the 

15 ECG complex, [2] the midpoint between the start of the ECG complex and the 
beginning of the QRS complex, [3] the beginning of the QRS complex, [4] the 
midpoint between the beginning of the QRS complex and the end of the QRS 
complex, [5] the end of the QRS complex, [6] the quarter point between the end of 
the QRS complex and the end of the ECG complex, [7] the midpoint between the 

20 end of the QRS complex and the end of the ECG complex, [8] the three-quarter 
point between the end of the QRS complex and the end of the ECG complex, and 
[9] the end of the ECG complex. Other segmentation techniques, such as using 
fixed time intervals from the start of the ECG complex, will produce acceptable 
results. The segment definitions are inputs that determine coefficients 

25 component.. The determine coefficients component generates transformation 
coefficients to be used in the synthesize of Lead Y. Other segmentation 
techniques, such as using time intervals generally corresponding to the start of the 
ECG complex, will produce acceptable results. In a preferred embodiment, the 
coefficients for a linear solution are generated using least squares analysis. 

30 Although the solution to the linear equation produces acceptable results, 
polynomial solutions also produce acceptable results. In addition, other curve- 
fitting solution may produce acceptable results. In an preferred embodiment, the 
outputs of determine coefficients (213) are the linear regression coefficients. Bo, 
Bj, Bo, and Bj. 

35 Figure 3 is a block diagram showing the components of Phase 2 in a 

preferred embodiment. The components illustrated are signal average (310), 
remove pacemaker spike (311), align ECG complex (312), and the generate 
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synthesized Y (313). In the signal average component, the system generates an 
averaged signal for each of Leads A, B, and C In a preferred embodiment, the 
system generates a running average of the last ten ECG complexes received. 
Various signal analysis techniques as described above for the signal average 
5 component (210) can be used. The signal average data is used as input to the 
remove pacemaker spike component and the align ECG complex component. 
The remove pacemaker spike component performs essentially the same function 
as described above for the remove pacemaker spike (211). The ECG complex 
component ensures that the ECG complex is aligned with the data processed in 
10 Phase 1. This alignment ensures that the coefficients will be applied to the proper 
segments. In the generate synthesized Y component, the system generates a 
synthesized Lead Y (Lead Y') based on the aligned signals for Leads A, B, and C, 
the coefficients, and the segment definition. In a preferred embodiment, a 
synthesized point of data is generated by the following equation- 

15 

Y/' = Bo + Bj » A{ + B 2 ' B/ + B 3 * Q 

where i indicates the /th data point in the ECG complex, Bj indicates the 
coefficients for the segment that contains the ith data point, and A{, Bj, and Q 
20 indicate the signal averaged data for the ith data point for leads A, B, and C 

In a preferred embodiment, the present invention is implemented in 
a computer based system. However, the present invention can be implemented 
using discrete logic in whole or in pan. Figure 9 is a flow diagram of a computer 
routine that generates the transformation coefficients. The input to this routine is 

25 an array a>ntaining the digitized data for leads I, n, and V 2 and for the lead to be 
synthesized. In a preferred embodiment, the array contains 300 entries for each 
lead. The data is sampled at a rate of 250Hz. In a preferred embodiment, the 
data is gathered over ten R-R cycle periods. For each lead, the ECG's are aligned 
and a signal analysis technique is employed to ininimize noise on the signal. In a 

30 preferred embodiment, the median point for each of the ten R-R cycles is used as 
a substitute for the average. Alternatively, the average value of the ten data 
points can be used. The output of this routine is the coefficients for the linear 
transformation. The signal averaged data is stored in array <?cg[300][4], which 
contains data for the base leads, and arrays [300], which contains data for the lead 

35 to be synthesized (see blocks 10A12 and 10H13). Columns 1, 2, and 3 of array ecg 
contain the base lead data and column 0 contains all l's and is used to calculate a 
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dc-offset (y-intercept of the line defined by the array and each of the two leads) 
for each segment 

In block 910, the system invokes subroutine DetermineJPacemaker, 
which analyzes the leads to determine whether the data was collected from a 
5 patient with a pacemaker. Subroutine DetermineJPacemaker returns a flag that 
indicates whether a pacemaker is present and it returns an index in array ecg 
specifying the highest point of the pacemaker spike. In block 911, if the 
pacemaker flag is true, then the system continues at block 912 to call subroutine 
Remove_Pacemaker_Spike, else the system continues at block 913. The 

10 techniques for detecting and removing pacemaker spikes are well-known. 

In block 913, the system calls subroutine SegmentJBCG. 
Subroutine Segment JECG analyzes the data in array ecg to determine the 
segments, and returns the start and stop indices into array ecg, which define each 
segment boundary, and the number of segments in variable segment j:ounl. The 

15 start and stop indices are returned in an array. The techniques for detecting the 
boundaries for the segments as . described above are well known. In blocks 914 
through 917, the system executes a loop that calculates the coefficients for each 
segment When the loop completes, the coefficients are stored in array B. Array 
5 is an array with four columns (one for each base lead and one for the dc-offset) 

20 and a number of rows equal to the number of segments. In block 914, the system 
initializes the loop control variable i to 0. In block 915, the system calls subroutine 
Regression, which determines the coefficients for the specified segment and stores 
the coefficients in array B. Figures 10A through 101 comprise a flow diagram for 
the regression subroutine. In block 916, the system increments the loop variable L 

25 In block 917, if the loop variable i equals se&nent jount, then coefficients for all 
the segments have been generated and the routine is done, else the system 
continues at block 915 to generate coefficients for the next segment 

Figures 10A through 101 comprise a flow diagram of subroutine 
Regression. The input parameters to subroutine Regression are array ecg (which 

30 contains the data for leads I, H, and and IS m column 0), array y (which 
contains the data for the lead to be synthesized), and variables start and stop 
(which are indices into the data arrays delimiting the segment boundaries). The 
output of subroutine Regression is array J5, which contains the coefficients for the 
linear transformation for the lead to be synthesized. The coefficients for the 

35 linear transformation are defined by the following equation: 

B = (ecg t *eegy l *(ecg L *y) 
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where ecg 1 is the transpose of matrix ecg and A indicates matrix inversion. 

In Figure 10A. the system calculates the transpose of ecg between 
start and stop and stores the transpose in array ecg t. Array ecg t is a 4-by-300 
5 array. In block 10A10, the system initializes loop variable i to 0. in block 10A11 
the system initializes inner loop variable / to equal parameter stem, which is the' 
index to the start of the segment in ecg. In block 10A12, the system set ecgj[i]{j] 
equal to ecgU][i]. In block 10A13, the system increments loop variable/. In block 
10A14, if loop variable; is greater than parameter stop (which indicates the end of 

10 the segment), then the system continues at block 10A15, else the system loops to 
block 10A12 to continue with the transposition of the matrix. In block 10A15, the 
system increments loop variable i. In block 10A16, if loop variable i equals 4, ihen 
the matrix transposition is complete and the system continues at block 10B10 of 
Figure 10B, else the system loops to block 10A11 to continue with the 

15 transposition of the matrix. 

In Figure 10B, the system multiplies ecg t times ecg and stores the 
result in array ecgj_ecg, which is a 4-by-4 array. In block 10B10, the system 
initializes loop variable / to 0. Loop variable i controls the looping through each 
of the four rows in eegj. In block 10B11, the system sets loop variable k equal to 
20 variable i. Loop variable k controls the looping through each of the four columns 
of ecg. Loop variable k is initialized to i rather than 0 because the product of a 
matrix and its transpose is a matrix that is symmetric about the diagonal. In block 
10B12, the system sets variable sum equal to 0. In block 10B13, the svstem sets 
loop variable ; equal to parameter start. Loop variable ; controls the looping 
through the segment data. In block 10B14, the system increments variable sum by 
the product of ecgjm times ecg{j][k]. In block 10B15, the system increments 
loop variable j. In block 10B16, if loop variable / is greater than parameter stop 
then all the data for the segment has been processed and the system continues ai 
block 10B17, else the system loops to block 10B14. In blocks 10B17 and 10B18, 
the system sets the symmetrical entries in the resulting matrix. In block 10B17, the 
system set ecgj_ecg[[][k] equal to variable sum. In block 10B18, the system sets 
ecgj_ecg[k)[i] equal to variable sum. In block 10B19, the system increments loop 
variable k. In block 10B20, if loop variable k is greater than 3, then each of 
columns of array ecg has been processed and the system continues at block 10B21 
else the svstem loops to block 10B12. In block 10B21, the system increments loop 
variable i. In block 10B22, if loop variable / is greater than 3, then all the rows of 
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array ecgj have been processed and the system continues at block 10C10 of 
Figure IOC, else the system loops to block 10B11. 

Figures 10C through 10F are a flow diagram for the inversion of 
matrix ecgj^ecg. In a preferred embodiment, the system uses the Gauss-Jordan 
5 inversion technique. Array inv_ecgj_ecg is an 4-by-8 array. This matrix is used 
for calculating the inverted matrix and for storing the inverted matrix in columns 0 
through 3. In Figure 10C, the system initializes matrix invjcgj_ecg for the matrix 
inversion. Columns 0 through 3 are initialized to contain array ecg t ecg and 
columns 4 through 7 are initialized to contain the identity matrix, that is, to 
10 containers in the diagonal and 0*s elsewhere. In blocks 10C10 through 10C16, the 
system copies array ecgj_ecg into columns 0 through 3 of array inv_ecg t_ecg. In 
block 10C10 t the system initializes loop variable i to 0. In block 10C11, the system 
initializes loop variable / to 0. In block 10C12, the system sets invj>cgjj>cg[i][f] to 
ecgj_ecg[i]\f]. In block 10C13, the system increments loop variable /. In block 
15 10C14, if loop variable j equals 4, then all the columns for the row specified by 
variable i have been processed and the system continues at the next row in block 
10C15, else the system loops to block 10C12. In block 10C15, the system 
increments loop variable /. In block 10C16, if loop variable i equals 4, then all the 
rows have been processed and the system continues at block 10C17, else the 
20 system loops to process the next row at block 10C11. 

In blocks 10C17 through 10C25, the system sets columns 4 through 7 
of mvjcgj_ecg equal to the identity matrix. In block 10C17, the system initializes 
loop variable / to 0. In block 10C18, the system initializes loop variable ; to 4. In 
block 1QC19, if loop variable i plus 4 equals loop variable /, then the loop variables 
25 index the diagonal of the matrix and the system continues at block 10C21, else the 
system continues at block 10C20. In block 10C20, the system sets the non- 
diagonal element inv_ecgj_ecg[i\\j\ equal to 0. In block 10C21, the system sets the 
diagonal element invjcgjjcg[i\\j] equal to 1. In block 10C22, the system 
increments loop variable In block 10C23, if loop variable / equals 8, then all the 
30 columns for the specified row have been set and the system processes the next row 
at block 10C24, else the system loops to block 10C19. In block 10C24, the system 
increments loop variable /. In block 10C25, if loop variable / equals 4, then all the 
rows have been processed and the system continues at block 10D10 in Figure 10D t 
else the system loops to block 10C18. 
35 Figure 10D is a flow diagram of a portion of the matrix inversion 

employing partial pivoting that ensures that for each of columns 0 through 3 the 
diagonal value is greater than the values below the diagonal in array inv ecg t ecg. 
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The system swaps the rows of the array to ensure this. In block 10D10. the system 
initializes loop variable; to 0. In block 10D11, the system sets variable pivot equal 
to mvjcgjjcgm. In block 10D12, the system sets variable largest equal to 
variable/. In block 10D13, the system initializes variable / to equal variable; plus 
5 1. In block 10D14, if the absolute value of inv ecg t ecg[i\\j] is greater than the 
absolute value of variable pivot, then the element indexed by variables / and ; is 
greater than the other elements above it in the column but at or below the 
diagonal and the system continues at block 10D15, else the system continues at 
10D17. In blocks 10D15 and 10D16, the system resets variables pivot and largest. 
10 In block 10D15, the system sets variable pivot equal to invjcg t ecg[i\\j]. In block 
10D16, the system sets variable largest equal to variable!. fcTblock 10D17, the 
system increments variable i. In block 10D18, if loop variable i equals 4, then all 
the rows in the specified column have been processed and the system continues at 
block 10D19, else the system loops to block 10D14. In block 10D19, if variable 
15 largest equals variable ;, then the diagonal contains the largest element and the 
system continues at block 10D21, else the system continues at block 10D20. In 
block 10D20, the system swaps the elements in the rows indexed by variables; and 
largest. In block 10D21, the system increments loop variable ;. In block 10D22, if 
loop variable ; equals 4, then columns 0 through 3 have been processed and the 
20 system continues at block 10E10 in Figure 10E, else the system loops to block 

In Figure 10E, the system generates the inverted matrix in columns 
4 through 7 of array invjcgjjcg. In block 10E10. the system initializes loop 
variable / to 0, which controls the looping through the rows. In block 10E11 the 
25 system initializes loop variable ; equal to variable i plus 1. Loop variable ; 
controls looping through row numbers greater than that specified by variable L In 
block 10E12, the system divides the element at inv ecg t ecg[i]{f] by 
inv_ecgj_ecg[i\[[\. In block 10E13, the system increments loop variably. In block 
10E14, if loop variable ; equals 8, then the dividing of the elements in the 
30 specified row is complete and the system continues at block 10E15, else the system 
loops to block 10E12. In block 10E15, the system sets inv ecg t_ecg[m equal to 
1. In block 10E16, the system initializes loop variable k to 0, which is an index to 
the rows. In block 10E17, if variable k equals variable /, then the system continues 
at block 10E23, else the system continues at block 10E18. In block 10E18, the 
35 system initializes variable ; equal to variable i plus 1. In block 10E19. the system 
subtracts invjcgj_ecg[k][i\ times invjcgjjcgflU] from invjcg t ecg[k)\j). In 
block 10E20, the system increments loop variable ;. In block ~10E21. if loop 
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variable / equals 8, then the elements in the specified row have been processed 
and the system continues at block 10E22, else the system loops to block 10E19. In 
block 10E22, the system sets invjcgj_ecg[k][[] equal to 0. In block 10E23, the 
system increments loop variable k. In block 10E24, if loop variable k equals 4 
5 then each row has been processed and the system continues at block 10E25, else 
the system loops to block 10E17. In block 10E25, the system increments variable 
i. In block 10E26, if loop variable i equals 4, then the matrix inversion is complete 
and the system continues at block 10F10 of Figure 10F, else the system loops to 
10E11. 

10 & 1 Figure 10F, the system moves the inverse matrix from columns 4 

through 7 of invjcgjjcg to columns 0 through 3. In block 10F10, the system 
initializes loop variable i to 0. In block 10F11, the system initializes loop variable/ 
to 0. In block 10F12, the system sets mvjcgj_ecg[i\\j\ to invjcgj_ecg[i\\j +4]. In 
block 10F13, the system increments loop variable /. In block 10F14, if loop 
15 variable / equals 4, then each element in the specified row has been copied and 
the system continues at block 10F15, else the system loops to block 10F12. In 
block 10F15, the system increments loop variable L In block 10F16, if loop 
variable / equals 4, then all the rows have been moved and the system continues at 
block 10G10 in Figure 10G, else the system loops to block 10F11. 
20 Figure 10G is flow diagram showing the initialization of arrays 

ec gJJ>W, which will contain the product of array ecgj times array y, and B, which 
will contain the coefficients. In block 10G10, the system initializes loop variable i 
to 0. In block 10G11, the system initializes ecgjj[i\ to 0. In block 10G1Z the 
system initializes B[segment-count][i] to 0, where segment count holds the index 
25 for the current segment In block 10G13, the system increments loop variable L 
In block 10G14, if loop variable / equals 4. then each element in the arrays has 
been initialized and the system continues at block 10H10 in Figure 10H, else the 
system loops to block 10G11. 

Figure 10H is a flow diagram showing the matrix multiplication of 
30 ecgj times y. The result is stored in array ecgjj, which is a 4-element array. In 
block 10H10, the system initializes loop variable ./ to 0. In block 10H11, the 
system initializes variable sum to 0. In block 10H12, the system initializes loop 
variable; by setting it equal to variable start. In block 10H13, the system adds 
ecgj[i\\j] times ytf] to the variable sum. In block 10H14, the system increments 
35 loop variable In block 10H15, if loop variable j is greater than variable stop, 
then the system continues at block 10H16, else the system loops to block 10H13. 
In block 10H16, the system sets ecgjj[{] equal to variable sum. In block 10H17, 
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the system increments variable i. In block 10H18, if loop variable i equals 4 then 
the system continues at block 10110 of Figure 101, else the system loops to block 
10H11. 

_ Figure 101 is a flow diagram showing the matrix multiplication of 

> wjcgjjcg times ecgjj. In block 10110, the system initializes loop variable i to 
a In block 10111, the system initializes variable sum to 0. In block 10IP the 
system initializes loop variable j to 0. In block 10113, the system "adds 
invjcgj_ecg[m times ecgjjrfl to variable sum. In block 10114, the svstem 
mcrements loop variable /. In block 10115, if loop variable / equals 4, then the 
10 system continues at block 10116, else the system loops to 10113. In block 10116 
the system sets B[segment-count][i] equal to variable sum. In block 10117 the' 
system mcrements loop variable L In block 10118. if loop variable / equals 4 'then 
calculation of the coefficients is complete and subroutine Regression returns^ else 
the system loops to block 10111. 

15 Fi § ure 11 » a flow diagram of routine Synthesize Lead. The input 

to routine Synthesize_Lead is the array of coefficients B, the segment boundaries 
and die signal averaged data from the three base leads which is stored in array 
fea*[3H300]. The routine generates the synthesized data for one ECG complex 
and stores the result in array* which is a 300 element array. In block 1110 the 
system initializes loop variable / to 0. In block llll, the system determines in 
which segment the element index by variable i is in. Variable segment is set to that 
segment. Variable segment is used as an index into array B to select the 
coefficients for the appropriate segment In block 1112, the svstem sets y[ij equal 

7^TT m ?IUS B[5e8mem M *"» b*Wm Plus B[segmem)[2) times 
-5 leadim Plus B{segmem][3) times lead{2}[i}. The addition of B{segmem][0l which 
corresponds to dc-offset, aligns the segments. Without adding in the dc-ofiset, the 
synthesized data has unacceptable intersegment gaps. In an alternate 
embodiment, the inter-segment gaps are effectively removed by a multivariate 
regression to mean 0. In block 1113, the system increments loop variable L In 
.0 block 1114, if loop variable i equals 300, then the generation of the synthesized 
lead is complete; else the system loops to block 1111. 

Although the present invention has been described in terms of a 
preferred embodiment, it is not intended that the invention be limited to these 
embodiments. Modifications within the spirit of the invention wfll be apparent to 
those skilled in the art For example, generation of the coefficients using 
polynomial equations will produce acceptable results. Also, it will be apparent to 
one skilled in the an that the use of segmentation can be applied to population- 
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based synthesis to produce improved results. It will also be apparent to one 
skilled in the art that the segmentation techniques of the present invention can be 
applied to data gathered from both orthogonal leads and non-orthogonal leads. It 
will also be apparent to one skilled in the art that non-standard leads can be 
5 synthesized by the present invention. It will also be apparent to one skilled in the 
an that determination of the coefficients can be made after a period of collecting 
data for the base leads. Thus, in the case of a medical emergency, the base lead 
data can be collected, then the coefficients determined and applied 
retrospectively, to analyze the ECG. It will also be apparent to one skilled in the 

10 art that the base leads data can be collected and the synthesis can be performed in 
a batch mode rather than in real-time. It will also be apparent to one skilled in 
the an that the 12 lead ECG can be synthesized from only three leads. Initially, 
data for leads I, II, Vy v 3t v 4» ^6 are collected. Then the transformation 
coefficients for leads V\, V3, V4, V$ 7 and Vg are determined using the leads I, H 

15 and V2 as base leads. The data for leads Vi, V3, V4, V5, and V$ can be collected 
simultaneously by using five electrodes or can be collected serially by using one 
electrode and moving it to the appropriate chest positions. In the synthesis mode, 
only data for the base leads is collected. The data for leads m, a VR, a VL, and a 
VF are calculated. The data for leads VI, V3, V4, V5, and V6 are synthesized 

20 based on the transformation coefficients. The scope of the present invention is 
defined by the claims that follow. 
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1. A method for synthesizing data for a given lead for an 
electrocardiogram for a patient, comprising the steps of: 
selecting a plurality of base leads; 

gathering a first set of electrocardiographic data from the patient for the 

ele~ I" f ,6ad "» - at leas on 

electrocardiographic complex; 

data; generatin e a tnmsformation based on the first set of electrocardiographic 

k , f?^ 8 a SeC ° nd 861 ° f ele ^«rdiographic data from the patient for 
complied " iMerVaI COrreSP ° nding t0 " 16351 ~ ^—graphic 

***** « T lyiR l tnmsformation to ^ond set of electrocardiographic 
data to effect the synthesis of data for the given lead. 

2. lie method of claim 1 wherein the step of generating a 
transformation includes the steps of: S 

lead., fl nH , f Wdhi ^ ° f eleOTO ^0graphic data for the base 

leads and for the given lead into segments; and 

generating a separate transformation for each segment; and 

wherein the step of applying the transformation includes the steps of: 

logically dividing the second set of electrocardiographic data for the base 
leads into segments; and 

applying the transformation for each segment to the data in the 
corresponding segment of the second set of electrocardiographic data to effect the 
synthesis of the data for the given lead. 

3. The method of claim 2 wherein the step of logically dividing the 
first set of electrocardiographic data into segments and the step of logically dividing the 
second set of electrocardiographic data into segments divide the electrocardiographic 
data into three segments defined (1) from the beginning of the electrocardiographic 
complex to the beginning of the QRS complex, (2) from the beginning of the QRS 
comp ex to the ending of the QRS complex, and (3) from the ending of the QRS 
complex to the ending of the electrocardiographic complex. 
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4. The method of claim 2 wherein the step of logically dividing the 
first set of electrocardiographic data into segments and the step of logically dividing the 
second set of electrocardiographic data into segments divide the electrocardiographic 
data into eight segments defined (1) from the beginning of the electrocardiographic 
complex to the midpoint between the beginning of the electrocardiographic complex 
and the beginning of the QRS complex, (2) from the midpoint between the beginning of 
the electrocardiographic complex and the beginning of the QRS complex to the 
beginning of the QRS complex, (3) from the beginning of the QRS complex to the 
midpoint between the beginning of the QRS complex and the ending of the QRS 
complex, (4) from the midpoint between the beginning of the QRS complex and the 
ending of the QRS complex to the ending of the QRS complex, (5) from the ending of 
the QRS complex to the point one-fourth the way between the ending of the QRS 
complex and the ending of the electrocardiographic complex, (6) from the point one- 
fourth the way between the ending of the QRS complex and the ending of the 
electrocardiographic complex to the midpoint between the ending of the QRS complex 
and the ending of the electrocardiographic complex, (7) from the midpoint between the 
ending of the QRS complex and the ending of the electrocardiographic complex to a 
point three-fourth the way between the ending of the QRS complex and the ending of 
the electrocardiographic complex, and (8) from a point three-fourth the way between 
the ending of the QRS complex and the ending of the electrocardiographic complex to 
the ending of the electrocardiographic complex. 

5. The method of claim 2 wherein the step of logically dividing the 
first set of electrocardiographic data into segments and the step of logically dividing the 
second set of electrocardiographic data into segments includes the steps of: determining 
time intervals for typical electrocardiographic data for an electrocardiographic complex 
corresponding to preselected portions of the electrocardiographic complex, and 
logically dividing the first and second sets of data based on the determined time 
intervals. 

6. The method of claim 2 wherein the step of logically dividing the 
first set of electrocardiographic data into segments and the step of logically dividing the 
second set of electrocardiographic data into segments, each set of electrocardiographic 
data having a QRS compex and an electrocardiographic complex, divides the data into 
segments such that the data in the QRS complex is in a different segment than the other 
data comprising the electrocardiographic complex. 
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7. The method of claim 3 wherein the step of applying the 
transformation includes the step of adding a dc-offset to the synthesized data to effect 
tne removal inter-segment gaps. 

8. The method of claim 4 wherein the step of applying the 
transformation includes the step of adding a dc-offset to the synthesized data to effect 
the removal inter-segment gaps. 

9. The method of claim 5 wherein the step of applying the 
transformation includes the step of adding a dc-offset to the synthesized data to effect 
the removal inter-segment gaps. 

10. The method claim 3 where in the step of applying the 
transformation includes the step of adjusting the synthesized data so that the mean of 
the data is equal to zero. 

11. The method claim 4 where in the step of applying the 
transformation includes the step of adjusting the synthesized data so that the mean of 
the data is equal to zero. 

12. The method claim 5 where in the step of applying the 
transformation includes the step of adjusting the synthesized data so that the mean of 
tne data is equal to zero. 

v, , , , 13 ' ^ method of 3 wherein the base leads are a subset of the 
12 standard electrocardiographic leads. 

it a , , K ™ e method of da ™ 4 wherein the base leads are a subset of the 
12 standard electrocardiographic leads. 

^ „• ^ , ^ meth ° d ° f Cl3im 12 wherein the »*« ^ the standard 
electrocardiographic leads includes leads I, II, and V?. 

16. The method of claim 13 wherein the subset of the standard 
electrocardiographic leads includes leads I, n, and V 2 . 
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17. The methocf of claims 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 or 

16 wherein 

the step of generating a transformation includes the step of determining 
transformation coefficients by applying the equation: 

B^(ec^ x ecg)' 1 *(ec^ M y) 

where B represents a vector (B i ... B n ) such that B{ is the transformation coefficient for 
base lead i, ecg represents a matrix containing the data from the first set of 
electrocardiographic data for the base leads, ecgf represents the transpose of matrix ecg, 
"1 represents matrix inversion, and y is a vector containing the data from the first set of 
electrocardiographic data for the lead to be synthesized; and 

the step of applying the transformation includes the step of determining 
the synthesized data by applying the equation: 

y'=5/*L!+ - +S rt *L a 

where y* represents a vector containing the synthesized data, £/ represents the 
transformation coefficient for base lead i, and Lj represents a vector containing the 
second set of electrocardiographic data for base lead i. 

18. The method of claims 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 or 

16 wherein 

the step of generating a transformation includes the step of determining 
transformation coefficients by applying the equation: 

B = (ecg/ * ecgjY 1 * (ecg/ ■ y) 

where B represents a matrix such that J5/j is the transformation coefficient for base lead 
i for the ;th order, ecgj represents a matrix containing a polynomial expression of the 
data from the first set of electrocardiographic data for the base leads, ecg/ represents 
the transpose of matrix ecgj, _1 represents matrix inversion, andy is a matrix containing 
the data from the first set of electrocardiographic data for the lead to be synthesized; 
and 



WO 92/02171 



PCT/US91/05562 



19 



the step of applying the transformation includes the step of determining 
the synthesized data by applying the equation: " ewnniniqg 
m n 

j = i / = i 

where y* represents a vector containing the synthesized data, B if represents the 
^sformatton coefficient for base lead i for the /th order, U represents a L 0 
containing me second set of electrocardiographic data for base lead In represents! 
number of base leads, and m represents the order of the polynomial express^ 

e l e «rn^,rH- * A meth ° d f ° r ****** *a ^ a given lead for an 
electrocardiogram for a patient, comprising the steps of: 

selecting a plurality of base leads; 

gathering a first set of electrocardiographic data from more than one 

01 — — — - - - - 

generatiDg a separate transformation for each segment based on the first 
set of electrocardiographic data; 

the h,< , , 8a ? erin8 " 5eC0Dd Mt ° f ^^^graphic data from the patient for 
debase ,=ads for an mttr val corresponding to a, ieast one tomographic 

ieads into J£?J*** "* ~ ° f - for rhe hase 

applying the transformation for each segment to the data in the 
correspondmg segment of the second set of electrocardiographic data to effect the 
synthesis of data for the given lead. 

™„ t , 2< \- 71,8 method 1' wherein the step of logically dividing first 

set of electiocardtographic dam into segmenrs and fit. step of logically dividing *e 
*cond se, of elecrrocardiogrophic data into segments divide the elecmLdio^ 
tot tnto three segments, each ser of electrocardiographic data having a beginn^ o 
tomographic complex, a beginning of the QRS comple, I end^gTthJ 

denned m ^ "t 1 " el «™»<«c complex, the three segments 

defined (1) from the begtnning of the electrocardiographic complex to the begJung of 



WO 92/02171 



PCT/US91/05562 



20 

the QRS complex, (2) from the beginning of the QRS complex to the ending of the 
QRS complex, and (3) from the ending of the QRS complex to the ending of the 
electrocardiographic complex. 

21. The method of claim 19 wherein the step of logically dividing the 
first set of electrocardiographic data into segments and the step of logically dividing the 
second set of electrocardiographic data into segments divide the electrocardiographic 
data into eight segments, each set of electrocardiographic data having a beginning of 
the electrocardiographic complex, a beginning of the QRS complex, an ending of the 
QRS complex, and an ending of the electrocardiographic complex, the eight segments 
defined (1) from the beginning of the electrocardiographic complex to the midpoint 
between the beginning of the electrocardiographic complex and the beginning of the 
QRS complex, (2) from the midpoint between the beginning of the electrocardiographic 
complex and the beginning of the QRS complex to the beginning of the QRS complex, 
(3) from the beginning of the QRS complex to the midpoint between the beginning of 
the QRS complex and the ending of the QRS complex, (4) from the midpoint between 
the beginning of the QRS complex and the ending of the QRS complex to the ending of 
the QRS complex, (5) from the ending of the QRS complex to the point one-fourth the 
way between the ending of the QRS complex and the ending of the 
electrocardiographic complex, (6) from the point one-fourth the way between the 
ending of the QRS complex and the ending of the electrocardiographic complex to the 
midpoint between the ending of the QRS complex and the ending of the 
electrocardiographic complex, (7) from the midpoint between the ending of the QRS 
complex and the ending of the electrocardiographic complex to a point three-fourth the 
way between the ending of the QRS complex and the ending of the 
electrocardiographic complex, and (8) from a point three-fourth the way between the 
ending of the QRS complex and the ending of the electrocardiographic complex to the 
ending of the electrocardiographic complex. 

22. The method of claim 19 wherein the step of logically dividing the 
first set of electrocardiographic data into segments and the step of logically dividing the 
second set of electrocardiographic data into segments divide the data into segments of 
predetermine lengths corresponding to predetermined time intervals. 

23. An apparatus for synthesizing the electrocardiographic data 
received from an electrode placed on a patient's body surface comprising: 
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a first, second, third, and fourth electrode attachable to the body surface, 
the electrodes responsive to the electrical activity of the heart; 

an analog-to-digital convener operatively connected to the electrodes to 
convert the electrical signal of the electrodes to digital data; 

data storage means operatively connected to the converter for storing 
digital data corresponding to an electrocardiographic complex for each electrode; 

calculating means for generating transformation data based on the stored 
digital data for the first, second, and third electrodes to the fourth electrode; and 

synthesizing means for applying the transformation data to digital data 
from the first, second, and third electrodes to synthesize data for the fourth electrode. 

24. The apparatus of claim 23 additionally comprising segmenting 
means for dividing the stored digital data into segments. 

25. The apparatus of claim 24 wherein the segmenting means divides 
the stored digital data into three segments, the stored digital data having a begmning of 
the electrocardiographic complex, a begmning of the QRS complex, an ending of the 
QRS complex, and an ending of the electrocardiographic complex, the three segments 
defined (1) from the beginning of the electrocardiographic complex to the begmning of 
the QRS complex, (2) from the beginning of the QRS complex to the ending of the 
QRS complex, and (3) from the ending of the QRS complex to the ending of the 
electrocardiographic complex. 



26. The apparatus of claim 24 wherein the segmenting means divides 
the stored digital data into eight segments, the stored digital data having a beginning of 
the electrocardiographic complex, a begmning of the QRS complex, an ending of the 
QRS complex, and an ending of the electrocardiographic complex, the eight segments 
defined (1) from the beginning of the electrocardiographic complex to the midpoint 
between the beginning of the electrocardiographic complex and the beginning of the 
QRS complex, (2) from the midpoint between the beginning of the electrocardiographic 
complex and the beginning of the QRS complex to the beginning of the QRS complex, 
(3) from the beginning of the QRS complex to the midpoint between the beginning of 
the QRS complex and the ending of the QRS complex, (4) from the midpoint between 
the beginning of the QRS complex and the ending of the QRS complex to the ending of 
the QRS complex, (5) from the ending of the QRS complex to the point one-fourth the 
way between the ending of the QRS complex and the ending of the 
electrocardiographic complex, (6) from the point one-fourth the way between the 
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ending of the QRS complex and the ending of the electrocardiographic complex to the 
mtdpomt between the ending of the QRS complex and I eZ g of 2 
electrocardiographic complex, (7) from the midpoint between the ending of the QRS 
complex and the ending of the electrocardiographic complex to a point three-fourth the 
way between the ending of the QRS complex and the ending of the 
eiectrocardiographic complex, and (8) from a point three-fourth the way between the 
ending of the QRS complex and the ending of the electrocardiographic complex to the 
ending of the electrocardiographic complex. 

27. The apparatus of claim 24 wherein the synthesizing means 
comprises means for adding a dc-offeet to the synthesized data to remove inter-segment 
gaps. 

28. The apparatus of claim 24 wherein the synthesizing means 
comprises means for adjusting the synthesized data so that the means of the data is 
equal to zero. 

29. The apparatus of claims 23, 24, 25, 26, 27, or 28 wherein the 
transformation data includes transformation coefficients and wherein the calculating 
means comprises: 

means for deternuning transformation coefficients by applying the 

equation: 

B = {ec£ * ecg)' 1 • (ecg* *y) 

where B represents a vector (Bj ... B n ) such that 5/ is the transformation coefficient for 
base lead i, ecg represents a matrix containing the data from the first set of 
electrocardiographic data for the base leads, ecg 1 represents the transpose of matrix ecg, 
represents matrix inversion, and y is a vector containing the data from the first set of 
electrocardiographic data for the lead to be synthesized; and 
wherein the synthesizing means comprises: 

means for applying the transformation data according to the equation: 
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where y represents a vector containing the synthesized data. B t represents the 
transformation coefficient for base lead i, and Lj represents a vector containing the 
second set of electrocardiographic data for base lead i. 
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